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The anisotropy of angular distributions of emitted nucleons and light charged particles for the 
asymmetric reaction system, 40 Ar+ 19 'Au, at b=6fm and Ebeam=35, 50 and lOOMeV/u, are inves¬ 
tigated by using the Improved Quantum Molecular Dynamics model. The competition between the 
symmetry potential and Coulomb potential shows large impacts on the nucleons and light charged 
particles emission in projectile and target region. As a result of this competition, the angular distri¬ 
bution anisotropy of coalescence invariant Y(n)/Y(p) ratio at forward regions shows sensitivity to 
the stiffness of symmetry energy as well as the value of Y(n)/Y(p). This observable can be further 
checked against experimental data to understand the reaction mechanism and to extract information 
about the symmetry energy at subsaturation densities. 


INTRODUCTION 


The nuclear symmetry energy plays a crucial role for 
understanding not only nuclear structures and reactions 
but also astrophysical phenomena. However, theoreti¬ 
cal predictions for the density dependence of symmetry 
energy of nuclear matter show large uncertainties away 
from the normal density 0,0- Many efforts have been 
devoted to probe and constrain the symmetry energy at 
both subsaturation and suprastaturation densities by an¬ 
alyzing observables from nuclear reactions, nuclear struc¬ 
tures and neutron stars, such as: isospin diffusion|3 t 
0 , double neutron to proton ratio (DR(n/p)) 10|-17j. 
light charged particle flow HUSH, 7T-/7r+0i-Ef7 neu¬ 
tron skin 29H3ll|, Giant Dipole Resonance and Pygmy 


Dipole Resonance [32 h 34] and masses of Isospin Analog 
State 35], alpha decay[36(, nrass-radius relationship and 
gra vitational waves from merging neutron star binaries 
37, |38]. Up to now, a consensus on the constraints of 
symmetry energy at subsaturation density 0, (40| has 
been obtained. Nevertheless, differences among the con¬ 
straints of symmetry energy obtained from different mod¬ 
els or approaches still exist. Understanding these dif¬ 
ferences is a challenge in this field and stimulates nu¬ 
clear physicists to propose new probes that are sensitive 
enough to further discriminate the stiffness of symmetry 
energy and to perform new experiments for a compre¬ 
hensive understanding of the mechanism of neutron-rich 
heavy ion collisions. 


Some probes related to the slower reaction process for 
semi-peripheral collisions in asymmetric reaction systems 
near the Fermi energy have been investigated or mea¬ 
sured to understand the reaction mechanism of neutron- 
rich heavy ion collisions(HIC) and to extract the infor¬ 


mation about symmetry energy [41H45I ]. For example, in 
reference j4lj, the isotopic composition of intermediate 
mass fragments (IMFs) emitted at mid-rapidity in semi¬ 
peripheral collisions of 124 Sn+ 64 Ni was analyzed and it 
was found that the IMFs emitted in the early stage of the 
reaction show larger values of < N/Z > and stronger an¬ 
gular anisotropy. The linear density dependence of sym¬ 
metry energy has been obtained by comparing the data 
with SMF model. Calculations from transport modeljifl] 
also show that the dynamical emission of nucleons and 
light charged particles (LCPs) for asymmetric reaction 
systems, such as Y{n)/Y{jp) ratios and Y(t)/Y( 3 He) ra¬ 
tios as a function of impact parameters, are also sensi¬ 
tive to the density dependence of symmetry energy due 
to the isospin migration/diffusion mechanism in neck re¬ 
gion. Very recently, R. S. Wang et al. analyzed the 
energy spectra of light charged particles in coincidence 
with fission fragments in 40 Ar+ 197 Au at 35 MeV/u [3 ]. 
The triggering condition on fission events in experiment 
has a bias on the semi-peripheral collisions. By defin¬ 
ing the ratio of the isotopic yield (normalized to that 

of proton) at 80° and 158°, R iao {X) = ( [Y(xj/Y(p)t aa )’ 
where X represents all species with Z=1 and 2, it is 
found that Ri SO increases with N/Z of the species, in¬ 
dicating a relatively high neutron richness of the emit¬ 
ted particles at smaller angle. According to the studies 


47j, this result can be attributed to the hierarchical 


feature of the semi-peripheral collisions at Fermi ener¬ 
gies, saying that the heavier fragments are faster than 
the lighter fragments along the beam direction and the 
Coulomb potential is stronger near the projectile/target 
region compared to that at midrapidities. Thus, one can 
expect different Y(n)/Y(p) ratio at project/target region 
and midrapidities, naturally leading to a non-unity value 
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of Riso- The different form of symmetry energy may 
leads to different values of anisotropy of isospin contents, 
R iso {X) = R x {91)/Rx(92), in the simulations. Further¬ 
more, the non-unity value of Ri SO also suggests the reac¬ 
tion process is an nonequilibrium process. An investiga¬ 
tion on angular distribution of the N/Z of the light parti¬ 
cles may reveal the sequence of the particle emission with 
different isospin composition, which carries the informa¬ 
tion of the symmetry energy at subsaturation densities. 


In this work, we mainly discuss the dynamical emis¬ 
sion of nucleons and LCPs at early stages for asymmetric 
nuclear reaction system 40 Ar+ 19 'Au at semi-peripheral 
collisions by using the ImQMD05 code. In our simu¬ 
lations, the beam energy dependence of angular distri¬ 
bution of coalescence invariant neutron to proton ratios 
is also discussed from 35MeV/u to lOOMeV/u. Based 
on our calculations, we propose that the anisotropy of 
angular distribution of coalescence invariant neutron to 
proton ratios may be used to understand the competi¬ 
tion between the Coulomb and symmetry potential, and 
may be used for probing the isospin effects and form of 

I 


symmetry energy. 


BRIEF INTRODUCTION OF IMPROVED 
QUANTUM MOLECULAR DYNAMICS MODEL 


Within Improved quantum molecular dynamics model 
(ImQMD05 code) [IH, 48, 49|, nucleons are represented 
by Gaussian wavepackets and the mean fields acting on 
these wavepackets are derived from energy density func¬ 
tional with the potential energy U that includes the 
Skyrme potential energy with just the spin-orbit term 
omitted: 


U — Up + Umd + U cou i (1) 

where U p is the Skyrme potential energy, U m d is the mo¬ 
mentum dependent potential energy, and U cou i is the 
Coulomb energy. The nuclear contributions are repre¬ 
sented in a local form with 
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Here, 5=(p n — p P )/{p n + Pp) is the isospin asymmetry, 
and p n , p p are the neutron and proton densities, respec¬ 
tively. A symmetry potential energy density of the form 


Qf(p/p 0 )' yi S 2 p is used in the following calculations. The 
energy density associated with the mean-field momentum 
dependence is represented by 


Umd = ^ r [ d 3 Pid 3 P2fNi(pi)fN2{P2)l-57[ln{l + 5 x 10 4 (Ap) 2 )] 2 

P 0 AT 1 J 


(4) 


where are nucleon phase space distribution functions, 
Ap= \p\ — p* 2 1, the energy is in MeV and momenta are 
in MeV/c. The resulting interaction between wavepack¬ 
ets is described in Ref. [501]. In this work, a=-356 MeV, 
/3=303 MeV, 77 = 7 / 6 , g sur =19A7 MeVfrn 2 , g SU r,iso=- 
11.35 MeVfm 2 , ^=35.19 MeV, and g p , T =0 MeV. The 
fragments are recognized by the isospin dependent mini¬ 
mum spanning tree method [51 1. 


potential for peripheral collisions for 40 Ar+ 197 Au by us¬ 
ing the ImQMD model. The time evolution of density, 
symmetry potential and single particle potential felt by 
neutron and proton at overlapped region are analyzed. 
The competition can be clearly observed at projectile 
and target region which leads to the different values of 
anisotropy of angular distribution of light particles for 
different form of symmetry potential. 


RESULTS AND DISCUSSION 


In this section, we examine the competition between 
the different form of symmetry potential and Coulomb 













3 


Reaction dynamics and time evolution of Vg yrn 
V q at neck region 


and 


Figure 1 shows the density contour plots for 
40 Ar+ 19 ‘Au reaction at the beam energy of 35MeV/u 
and impact parameter b= 6 fm. In the calculations, the 
projectile and target start to touch at around 50fm/c, 
and then the overlapped region is formed and it reaches 
maximum compression at around 90fm/c. After about 
150fm/c, the overlapped region expands to lower density 
and ruptures into fragments. For semi-peripheral colli¬ 
sions we studied, the neck dynamics play important roles 
for the dynamical emission of LCPs[52j, and its yields and 
isospin contents for nucleons and LCPs are closely related 
to the symmetry potential in the participant region of 
the reaction. Thus, it is useful for us to understand how 
the symmetry potential and single particle potential for 
neutron and proton in the participant region evolve with 
time, because the single particle potential and symmetry 
potential will determine the yields and isospin contents 
of nucleons and LCPs emission. 



z (fm/c) 

FIG. 1. (Color online) Time evolution of the density con¬ 
tour plots in the reaction plane for 40 Ar+ 197 Au collisions at 
Ebeam =35 MeV/u and b=6fm from one typical event; the 
blue points represent the overlapped region of projectile and 
target. 

Figure 2 (a) shows the time evolution of the symmetry 
potential V£ ym = %-[(7 — l)it 7 5 2 ± 2 u^S], q = n,p , at 
the overlapped region which is presented as blue points 
in Figure 1. The point is defined as the equal distance 
from the surface of the projectile and the target, and a 
spherical region with two times of nucleon radius, i.e. 
R=1.7fm, is used to approximately describe the cen¬ 
ter of overlapped region in the calculations. The solid 
lines are the results from 7i=0.5 and dashed lines are for 
7i=2.0. The magnitude of the symmetry potential felt by 
neutron (V^ m ) for 7j=0.5 is obviously larger than that 
with 7 i= 2.0 except for the time at around 80fm/c where 
V s ™ m (7* = 0.5) < = 2.0). The reason is that the 

compressed density is close and up to normal density at 
around 80fm/c, and its isospin asymmetry comes close to 
0.2 due to isospin migration. Thus, the symmetry poten¬ 
tial felt by neutron (V s ™ m ) obtained with 7i=2.0 becomes 
larger than that for 7i=0.5 at p > 0.9po based on the for¬ 
mula of Vg ym [Ill HI] ■ However, the magnitude of symme¬ 
try potential felt by proton for 7 ^= 2.0 is still smaller than 
that for 7i=0.5 at p < l.lpo for S ~ 0.2, and it leads to 


the larger magnitude of symmetry potential felt by pro¬ 
ton for 7 j= 0.5 than that for 7 ^= 2.0 at different time for 
the beam energy of 35MeV/u. Since HIC observables ob¬ 
tained in the simulations are determined by the reaction 
process from early stages to late stages, one can expect 
that yield ratio of neutron to proton, i.e. Y(n)/Y(p), 
obtained with 7 j= 0.5 is greater than that with 7 ^= 2.0 
for neutron rich reaction system. Quantitatively, under¬ 
standing the Y(n)/Y(p) ratios requires the knowledge of 
single particle potential felt by neutron and proton in the 
reaction system. Figure 2 (b) and (c) show the the single 
particle potential for neutron and proton for 7 j= 0.5 and 
2.0, respectively. The solid lines are for protons where 
the Coulomb contribution is included, i.e., V p + Vcoui- 
The dashed lines are the nucleonic potential felt by the 
proton, i.e. V p , without Coulomb potential. The dotted 
lines are the single particle potential for neutrons, i.e. 
V n . Figure 2 (b) shows that V p + Vcoui is close to V n 
for 7 j= 0 . 5 . However, for 7i=2.0, V p + Vcoui is obviously 
larger than V n after 120fm/c because symmetry poten¬ 
tial becomes weak, especially at lower density. Thus, the 
Coulomb potential obviously moves the single particle 
potential felt by proton up to the strength by neutron. 
It may lead a higher yield for protons than neutrons, 
if the Coulomb repulsion is strong enough. Especially 
for asymmetric reaction systems at semi-peripheral colli¬ 
sions, one can expect to observe this competition between 
the Coulomb and symmetry potentials which influences 
the isospin contents for the emitted nucleons and LCPs, 
and hence, their angular distribution anisotropies show 
the sensitivity to the different forms of symmetry energy. 



t (fm/c) 


FIG. 2. (Color online) Time evolution of symmetry potential 
for neutron and proton obtained with 7 ^= 0 .5 (solid lines) and 
2.0 (dashed lines)(a); single particle potential for neutron and 
proton for 7 ; = 0.5 (b); and 7 i = 2.0 (c) at the overlapped 
region for 40 Ar-|- 197 Au at b= 6 fm and A(, eam =35MeV/u. 
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Angular distribution of the coalescence invariant 
yields for neutron and proton 


As is known, the absolute yields of LCPs are not 
very well reproduced due to the considerable deviation of 
binding energy of light clusters between the QMD model 
and experimental data 54J. Loosely bound clusters like 
d , t are over predicted, whereas strongly bound clusters 
like 4 He is underpredicted due to its lower binding energy 
in the QMD model. Therefore, it is hard to draw a firm 
conclusion by comparing the absolute number of these 
light fragments to data to extract the information about 
physical quantities, such as the equation of state. One 
method to eliminate the problem related to the absolute 
yield of LCPs is to introduce the coalescence invariant 
neutron and proton yields, and to construct the related 
observables, such as Coalescence Invariant Y(n)/Y(p) 
(Cl n/p ratio) ratio spectra 11, 12, [H^] . In this paper, 
we investigate the angular distributions of Cl n/p ratio 
by varying the stiffness of symmetry potential. The an¬ 
gular distribution of coalescence invariant (Cl) neutron 
and proton yields are constructed by adding the neutrons 
and protons in light particles to free neutrons and pro¬ 
tons at given 6 c , m , as follows: 


position does not obviously vary with time. As shown in 
Figure 3, the nucleons emitted at 9 cm . < 25° are mainly 
from the projectile while those at 9 cm . > 60° are mainly 
from the target. The emission of nucleons within the 
angular gate 20° < 9 c .m. < 60° corresponds to the com¬ 
pressed overlap region between the projectile and target. 
The anisotropy of angular distribution of Cl nucleons is 
due to the asymmetry of projectile and target in semi¬ 
peripheral collisions. Furthermore, one also observes that 
the yields of neutrons are obviously larger than that for 
protons after 150fm/c for 7j=0.5. But for 7i=2.0, the 
yields of protons from target are always higher than that 
for neutrons when 0 c . m . > 60°. The reason is that the 
strength of the symmetry potential is weak for 7 j= 2.0 
at low density regions where many nucleons are emitted, 
and the target-like fragments provide stronger Coulomb 
repulsive. Thus, the yields of protons become larger than 
neutrons for 7^=2.0. For the nucleons from the projectile, 
one also observe that yields of protons are greater than 
neutrons before 300fm/c, but it finally turns over. Af¬ 
ter 300fm/c, the yields of protons and neutrons become 
similar because of the weaker Coulomb potential in the 
projectile region. 


dM n>C i ^ Ar dY(N,Z) 

-Ta - = l^ N -—ra - W 

W't'c.m. jy g uvc.m. 

dM p , C i v dY(N,Z) 

dd c .m. j^ z d0 c . m . ’ 1 > 

the summation is over n, p, d, t, 3 He, 4 He and 6 He in this 
work, and Y(N,Z) is the yield of fragments with charge 
number Z and neutron number N. 

In order to understand the origin of angular distribu¬ 
tion of Cl yields for neutron and proton in asymmetric 
reaction systems 40 Ar+ 197 Au at 35MeV/u, we show the 
angular distribution of the yields for Cl neutrons and Cl 
protons emitted from the projectile (line with circles) and 
the target (line with squares) at 150fm/c, 200fm/c and 
350fm/c in Figure 3. Left panels are the results for 7i=0.5 
(Fig.3(a) and (c)), and right panels are for 2.0 (Fig.3(b) 
and (d)). The line with open (solid) symbols refers to the 
yields of Cl protons in Fig.3 (a) and (b) (Cl neutrons are 
in Fig.3(c) and (d)). The calculated results show that the 
yields of Cl nucleons from target have a wide angular dis¬ 
tribution and exhibit weak anisotropy in the angular dis¬ 
tributions. The peak of angular distribution appears be¬ 
low 9c.m. ~ 90° when the time is earlier than ~300fm/c, 
i.e., the nucleons and LCPs prefer to emit at forward re¬ 
gions for asymmetric reaction systems because the neck 
dynamics dominate the nucleons’ and LCPs’ emission at 
earlier stages. After then, the peak moves to the back¬ 
ward regions at late stages because the evaporation from 
target-like fragments becomes more and more important. 
The yields of Cl nucleons from the projectile obviously 
peak around 9 c m . ~ 25° and the corresponding angular 



FIG. 3. (Color online) The angular distribution of the yields 
for Cl neutrons ((a) and (b)), and Cl protons ((c) and (d)) 
emitted from projectile (line with circle symbols, left curves in 
each panel) or target (line with square symbols, right curves 
in each panel) at 150fm/c, 200fm/c and 350fm/c. Left panels 
are for 7 ; = 0.5 and right panels are for 2.0. 

Figure 4(a) shows the angular distributions of yields 
of Cl neutrons (solid symbols) and protons (open sym¬ 
bols) for 7j=0.5 (red squares) and ji=2.0 (blue cir¬ 
cles) at the stop time of simulations, t=400fm/c, for 
Ebeam= 35MeV/u and b=6fm. In general, the yields of 
nucleons obtained with 7i=0.5 are greater than that for 
7 i= 2.0 due to its stronger symmetry energy at subsatu¬ 
ration densities. For the isospin contents, the yields of 
neutrons are always greater than the yields of protons at 
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whole angular range for 7 j= 0 . 5 . Similar to the results at 
35MeV/u, the yields of neutrons are greater than protons 
for the beam energy at 50 (Figure 4 (b)) and lOOMeV/u 
(Figure 4 (c)). But for the 7 j= 2 . 0 , a different behavior 
appears in the yields of neutrons and protons. At the 
beam energy of 35MeV/u, the yields of neutrons become 
less a little bit than the yields of protons at 9 c .m. > 60° 
for 7 j= 2 . 0 . It is the result of the competition between the 
Coulomb potential and symmetry potential for asymmet¬ 
ric reaction system at semi-peripheral collisions as men¬ 
tioned before. At the beam energy of lOOMeV/u, the 
reaction is more violent and a larger part of the collid¬ 
ing system is dissociated into gases (nucleons and light 
particles), and thus more neutrons originally bounded in 
the system are released and outnumber the protons for 
the whole angular region, as shown in Figure 4(c). 



FIG. 4. (Color online) Angular distribution of the yields for 
Cl neutrons (solid symbols) and Cl protons (open symbols) 
emitted at t=400fm/c obtained with 7 i = 0.5 (red symbols) 
and 2.0 (blue symbols). The left panel is the result at the 
beam energy of 35MeV/u, middle panel is for 50MeV/u and 
right panel is the result for lOOMeV/u. The impact parameter 
is 6 fm. 


Correlation between the symmetry energy and 
angular distribution of Cl n/p ratios 

Consequently, we construct the angular distribution of 
Cl n/p ratios, i.e. d gg n '^ 1 / d dg P ' CI ratios, to explore the 
information of reaction mechanism and symmetry energy. 
As shown in Figure 5, the Cl n/p ratios obtained with 
7 i= 0.5 are greater than that with 7 ,=2.0 for the beam en¬ 
ergy we studied, because the symmetry energy for 7 j= 0.5 
is stronger than 7 ^= 2 .0 at subsaturation density. The 
sensitivity of Cl n/p ratios to 7 i becomes large at the 
target region where the charge number of fragment is 
large, and it is consistent with the results in [56| where 
the influence of symmetry potential and in-medium NN 
cross section on the free nucleon n/p ratio were studied 


for the central collision of 100 Zn-|- 40 Ca at 200MeV/u. 

In addition to the values of Cl n/p ratio, analyzing 
the d c .m. dependence on the Cl n/p ratio will provide 
more information for revealing the reaction mechanism 
on the competition between the Coulomb and symmetry 
potential. For £/ eaTO =35MeV/u, the Cl n/p ratios show 
a different 0 c . m . dependence for the different forms of 
symmetry potential. In case of 7 ^= 0 .5, the Cl n/p ratios 
slightly increase as a function of 0 c . m .. It is the result of 
the isospin asymmetry changing from d pr - O j= 0.10 at the 
projectile region to 5tar= 0.19 at the target region since 
the single particle potentials of neutrons and protons are 
close, as shown in Figure 2 (b). However, in case of 
7 i= 2 . 0 , the Cl n/p ratios obviously decrease with 0 c . m ., 
and the Cl n/p ratios at backward are smaller than that 
at forward regions. Furthermore, we also analyze the an¬ 
gular distributions of the Cl n/p ratios at beam energy 
for 50 and lOOMeV/u. The Cl n/p ratios increase as 
the beam energy increases except for the Cl n/p at for¬ 
ward region and 7 i= 0 . 5 . The increasing behavior of the 
Cl n/p ratios as the beam energy increases are from the 
enhancement of the symmetry potential at larger com¬ 
pressed density which can be achieved with higher beam 
energies. However, Cl n/p ratios at forward regions for 
7 i= 0.5 show opposite behaviors, i.e. the Cl n/p ratios 
decrease with the beam energy increasing. It is the ef¬ 
fects of the finite size of the reaction systems. The final 
values of the Cl n/p ratios from heavy ion collisions not 
only depend on the symmetry potential, but also depend 
on the yield of nucleons for finite reaction system. Since 
there are too many nucleons and LCPs emitted from the 
projectile for 7 i= 0.5 at Sb eam =100MeV/u, the Cl n/p 
ratios at forward regions tend to decrease to the initial 
N/Z value of projectile, N/Z pro j= 1.22, and is obviously 
smaller than that obtained at 35MeV/u. 

To check the smearing effects of impact parameters on 
the sensitivity of this observable to 7 j, calculations for 
b=4 and 8 fm have also been performed with ImQMD. 
The calculations show that the Cl n/p ratio insignifi¬ 
cantly depends on the impact parameters from b=4 to 
8 fm at the beam energy of 35, 50 and 100 MeV/u, be¬ 
cause the value of Cl n/p ratio reflects the symmetry 
energy information of reaction system which is not var¬ 
ied dramatically from b=4 to 8 fm for 40 Ar+ 19 ' Au. For 
its insensitivity to the changing of the impact parameter, 
the Cl n/p ratio for mid-peripheral collisions shall be a 
valid probe to the symmetry energy. 

Focusing on the anisotropy of angular distributions 
of the Cl n/p ratios and its sensitivity to 7 ,; (or the 
slope of symmetry energy, L = 3po^§^U;)> we an¬ 
alyzed R iso 1 = i? G7n/p (20 o )/i? c . / „/ p (90 o ) and R iso2 = 
i? G / rl / p (90 o )/i? G / rl / p (160 o ), which makes the experimen¬ 
tal discrimination easier, for b= 6 fm at different beam en¬ 
ergies. The calculations show R iso 1 increases by ~20% 
when 7 i varying from 0.5 to 1.5 for £'h eQm =35MeV/u. 
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FIG. 5. (Color online) Left panel shows the calculated results 
for Cl n/p ratios as a function of # c . m . for b=4, 6 and 8 fm 
at i?i,eam=35MeV/u. The solid symbols are for 7 ; = 0.5 and 
open symbols are for 7 ; = 2.0. Middle panel is the results for 
50MeV/u and right panel is for lOOMeV/u. 

For Riso 2 ; which is constructed from the backward re¬ 
gions, its values weakly depend on the symmetry energy 
for the beam energies we studied. Similar behaviors have 
also been found for 50 and lOOMeV/u. These calcula¬ 
tions illustrate that Ri SO \ is strongly correlated to the 
stiffness of symmetry energy, and it would be helpful for 
us to extract the information of symmetry energy and to 
understand the reaction mechanism as well as the spectra 
of Cl n/p ratios if more complete data are measured. 

pMeV; L(MeV) 



FIG. 6 . (Color online) The calculated results of Ri so i 
(left panel) and Ri so 2 (right panel) as a function of 7 ; for 
40 Ar+ 197 Au at b= 6 fm and the beam energy is 35, 50 and 100 
MeV/u. 


SUMMARY 

In summary, the semi-peripheral heavy-ion collisions 
of 40 Ar-|- 197 Au at Eb eam =35, 50 and lOOMeV/u were 
studied by means of the improved quantum molecular 
dynamics model (ImQMD05). The competition between 


the symmetry potential and Coulomb potential is dis¬ 
cussed for the asymmetric reaction system. The calcula¬ 
tions show that the impacts of the competition between 
Coulomb and symmetry potential on the heavy ion col¬ 
lision observables, Cl n/p ratio, can be observed at the 
projectile and target regions. In case of 7 ;= 2 . 0 , the yields 
of protons are larger than neutrons in the target region 
for beam energy at 35MeV/u due to its stronger Coulomb 
potential and weaker symmetry potential. As the beam 
energy increasing up to lOOMeV/u, the reaction is more 
violent and a larger part of the colliding system is disso¬ 
ciated into gases (nucleons and light particles), and thus 
more neutrons originally bounded in the system are re¬ 
leased and outnumber the protons. For the 7 ^= 0 .5, the 
yields of neutrons are obviously greater than protons due 
to its stronger symmetry potential at subsaturation den¬ 
sity. The competition between the Coulomb and sym¬ 
metry potential causes the anisotropy of angular distri¬ 
bution of coalescence invariant n/p ratio, i?,: so i, which 
is accordingly sensitive to the stiffness of symmetry en¬ 
ergy as well as the values of Cl n/p ratio. For 7 ^= 2 .0, 
the Cl n/p ratios decrease with angle increasing, and 
Risoi = 1.14±0.01. For 7 i= 0 . 5 , the Cl n/p ratios slightly 
increase with angle because the single particle potential 
felt by proton is close to neutron, and R iso 1 = 0.94±0.01. 
The values of R iso 1 increase by about ~20% when 74 
varies from 0.5 (L=51MeV) to 1.5(L=104MeV). The sen¬ 
sitivity of Risoi to 7 i remains at all the beam energies we 
studied. It is of high interest to measure the isospin com¬ 
position of the light particles over a wide angular range 
in further experiments at Fermi energies to gain com¬ 
prehensive understanding on the reaction mechanism as 
well as on the the symmetry energy behavior at subsatu¬ 
ration densities, and the symmetry energy in warm dilute 
matter. 
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